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A quantum system consisting of two subsystems is separable if its den- 
sity matrix can be written as p = J^wkp'k ® P'ky where p'^ and p'^ 
are density matrices for the two subsytems, and the positive weights 
wk satisfy J^'f^K = 1- A necessary condition for separabihty is derived 
and is shown to be more sensitive than Bell's inequality for detecting 
quantum inseparability. Moreover, collective tests of Bell's inequality 
(namely, tests that involve several composite systems simultaneously) 
may sometimes lead to a violation of Bell's inequality, even if the latter 
is satisfied when each composite system is tested separately. 

1. INTRODUCTION 

From the early days of quantum mechanics, the question has often 
been raised whether an underlying "subquantum" theory, that would 
be deterministic or even stochastic, was viable. Such a theory would 
presumably involve additional "hidden" variables, and the statistical 



predictions of quantum theory would be reproduced by performing suit- 
able averages over these hidden variables. 

A fundamental theorem was proved by Bell [||] , who showed that if 
the constraint of locality was imposed on the hidden variables (namely, 
if the hidden variables of two distant quantum systems were themselves 
be separable into two distinct subsets), then there was an upper bound 
to the correlations of results of measurements that could be performed 
on the two distant systems. That upper bound, mathematically ex- 
pressed by Bell's inequality [0], is violated by some states in quantum 
mechanics, for example the singlet state of two spin-| particles. 

A variant of Bell's inequality, more general and more useful for 
experimental tests, was later derived by Clauser, Horne, Shimony, and 
Holt (CHSH) 0. It can be written 

I (AB) + {AB') + (A'B) - {A'B') \ < 2. (1) 

On the left hand side, A and A' are two operators that can be mea- 
sured by an observer, conventionally called Alice. These operators do 
not commute (so that Alice has to choose whether to measure A or A') 
and each one is normalized to unit norm (the norm of an operator is 
defined as the largest absolute value of any of its eigenvalues). Like- 
wise, B and B' are two normalized noncommuting operators, any one 
of which can be measured by another, distant observer (Bob). Note 
that each one of the expectation values in Eq. (|l]) can be calculated by 
means of quantum theory, if the quantum state is known, and is also 
experimentally observable, by repeating the measurements sufficiently 
many times, starting each time with identically prepared pairs of quan- 
tum systems. The validity of the CHSH inequality, for all combinations 
of measurements independently performed on both systems, is a nec- 
essary condition for the possible existence of a local hidden variable 
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(LHV) model for the results of these measurements. It is not in general 
a sufficient condition, as will be shown below. 

Note that, in order to test Bell's inequality, the two distant ob- 
servers independently measure subsytems of a composite quantum sys- 
tem, and then report their results to a common site where that informa- 
tion is analyzed 0. A related, but essentially different, issue is whether 
a composite quantum system can be prepared in a prescribed state by 
two distant observers who receive instructions from a common source. 
For this to be possible, the density matrix p has to be separable into a 
sum of direct products. 



where the positive weights wk satisfy J^wr = 1, and where p'j^ and p'^ 
are density matrices for the two subsystems. A separable system always 
satisfies Bell's inequality, but the converse is not necessarily true [4-7]. 
I shall derive below a simple algebraic test, which is a necessary condi- 
tion for the existence of the decomposition (|^). I shall then give some 
examples showing that this criterion is more restrictive than Bell's in- 
equality, or than the a-entropy inequality 0. 

2. SEPARABILITY OF DENSITY MATRICES 

The derivation of the separability condition is easiest when the 
density matrix elements are written explicitly, with all their indices . 
For example, Eq. (0) becomes 



K 

Latin indices refer to the first subsystem, Greek indices to the second 
one (the subsystems may have different dimensions). Note that this 
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equation can always be satisfied if we replace the quantum density 
matrices by classical Liouville functions (and the discrete indices are 
replaced by canonical variables, p and q). The reason is that the only 
constraint that a Liouville function has to satisfy is being non-negative. 
On the other hand, we want quantum density matrices to have non- 
negative eigenvalues , rather than non- negative elements, and the latter 
condition is more difficult to satisfy. 
Let us now define a new matrix, 

^mij,,nu — Pn^,mu- (4) 

The Latin indices of p have been transposed, but not the Greek ones. 
This is not a unitary transformation but, nevertheless, the a matrix is 
Hermitian. When Eq. @ is valid, we have 

a = J2wA{p'Af®PA- (5) 

A 

Since the transposed matrices {p'^)^ = {p'a)* non- negative matrices 
with unit trace, they can also be legitimate density matrices. It fol- 
lows that none of the eigenvalues of a is negative. This is a necessary 
condition for Eq. (0) to hold 0. 

Note that the eigenvalues of a are invariant under separate unitary 
transformations, U' and U", of the bases used by the two observers. In 
such a case, p transforms as 

p^{U'®U")p{U'(^U")\ (6) 

and we then have 

a^{U'^®U")a{U'^®U")\ (7) 

which also is a unitary transformation, leaving the eigenvalues of a 
invariant. 
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As an example, consider a pair of spin-| particles in an impure 
singlet state, consisting of a singlet fraction x and a random fraction 



(1 — x) Note that the "random fraction" (1 — x) also includes 

singlets, mixed in equal proportions with the three triplet components. 
We have 

Pmfj,,nu Smfi,nu ~l~ (1 ^mn /4, (8) 

where the density matrix for a pure singlet is given by 

'S'01,01 = 'S'10,10 = = — 'S'10,01 = \i (9) 

and all the other components of 5* vanish. (The indices and 1 refer to 
any two orthogonal states, such as "up" and "down.") A straightfor- 
ward calculation shows that a has three eigenvalues equal to (1 + x)/4, 
and the fourth eigenvalue is (1 — 3x)/4. This lowest eigenvalue is pos- 
itive if X < |, and the separability criterion is then fulfilled. This 
result may be compared with other criteria: Bell's inequality holds for 
X < l/-\/2, and the a-entropic inequality for x < l/-\/3. These 
are therefore much weaker tests for detecting inseparability than the 
condition that was derived here. 

In this particular case, it happens that this necessary condition is 
also a sufficient one. It is indeed known that if x < it is possible to 
write p as a mixture of unentangled product states ||TT| . This suggests 



that the necessary condition derived above (a has no negative eigen- 
value) might also be sufficient for any p. A proof of this conjecture was 



indeed recently obtained ]12[ for composite systems having dimensions 
2x2 and 2x3. However, for higher dimensions, the present necessary 
condition was shown not to be a sufficient one. 

As a second example, consider a mixed state consisting of a fraction 
X of the pure state a|01) -|- 6|10) (with |ap + |6p = 1), and fractions 
(1 — x)/2 of the pure states |00) and |11). We have 
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Poo,oo = Pii,ii = (1 - a;)/2, (10) 

Poi,oi = x\a\^, (11) 

Pio,io = xlb]'', (12) 

Poi,io = pIo,oi = xo-b*, (13) 

and the other elements of p vanish. It is easily seen that the a matrix 
has a negative determinant, and thus a negative eigenvalue, when 

X > {1 + 2\ab\)-\ (14) 

This is a lower limit than the one for a violation of Bell's inequality, 
which requires [[7| 

X > [1 + 2\ab\{V2 - (15) 

An even more striking example is the mixture of a singlet and a 
maximally polarized pair: 

Pmfi,nu = X Smfi,nu + (1 — x) 6mO ^nO 5^0 ^uO- (16) 

For any positive x, however small, this state is inseparable, because a 
has a negative eigenvalue (— x/2). On the other hand, the Horodecki 



criterion ||T3[ gives a very generous domain to the validity of Bell's in- 



equality: X < 0.8. 



3. COLLECTIVE TESTS FOR NONLOCALITY 

The weakness of Bell's inequality as a test for inseparability is 
due to the fact that the only use made of the density matrix p is for 
computing the probabilities of the various outcomes of tests that may 
be performed on the subsystems of a single composite system. On the 
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other hand, an experimental verification of that inequahty necessitates 
the use of many composite systems, all prepared in the same way. 
However, if many such systems are actually available, we may also test 
them collectively, for example two by two, or three by three, etc., rather 
than one by one. If we do that, we must use, instead of p (the density 
matrix of a single system), a new density matrix, which is p ® p, or 
p ® p ® Pi in a higher dimensional space. It will now be shown that 
there are some density matrices p that satisfy Bell's inequality, but for 
which p ® p, oi p ® p ® p, etc., violate that inequality . 

The example that will be discussed is that of the Werner states 
defined by Eq. Let us consider n Werner pairs. Each one of the 
two observers has n particles (one from each pair). They proceed as 
follows. First, they subject their n-particle systems to suitably chosen 
local unitary transformations, f/, for Alice, and V , for Bob. Then, they 
test whether each one of the particles labelled 2, 3, . . . , ra, has spin up 
(for simplicity, it is assumed that all the particles are distinguishable, 
and can be labelled unambiguously). Note that any other test that 
they can perform is unitarily equivalent to the one for spins up, as this 
involves only a redefinition of the matrices U and V . If any one of 
the 2{n — 1) particles tested by Alice and Bob shows spin down, the 
experiment is considered to have failed, and the two observers must 
start again with n new Werner pairs. 

A similar elimination of "bad" samples is also inherent to any ex- 
perimental procedure where a failure of one of the detectors to fire 
is handled by discarding the results registered by all the other detec- 
tors: only when all the detectors fire are their results included in the 
statistics. This obviously requires an exchange of classical information 
between the observers. (There is a controversy on whether a violation 
of Bell's inequality with postselected data O] is a valid test for non- 
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locality |I6[. I shall not discuss this issue here; I only examine whether 
or not Bell's inequality is violated by the postselected data.) 

The calculations shown below will refer to the case n = 3, for defi- 
niteness. The generalization to any other value of n is straightforward. 
Spinor indices, for a single spin-^ particle, will take the values (for 
the "up" component of spin) and 1 (for the "down" component). The 
16 components of the density matrix of a Werner pair, consisting of a 
singlet fraction x and a random fraction (1 — x), are, in the standard 
direct product basis: 

Pmn,st ^mn,st ~l~ (^ •^) ^ms ^nt /4, (17) 

where I am now using only Latin indices, contrary to what I did in 
Eq.(^; this is because Greek indices will be needed for another purpose, 
as will be seen soon. Thus, now, the indices m and s refer to Alice's 
particle, and n and t to Bob's particle. 

When there are three Werner pairs, their combined density matrix 
is a direct product p p' (g) p", or exphcitly, Pmn,st Pm>n' ,s't' Pm"n" ,s"t" ■ 
The result of the unitary transformations U and V is 

p®p' ®p" ^{U® V) {p®p® p") {U^ ® V^). (18) 

Explicitly, with all its indices, the U matrix satisfies the unitarity rela- 
tion 

mm' m" 

In order to avoid any possible ambiguity, Greek indices (whose values 
are also and 1) are now used to label spinor components after the 
unitary transformations. Note that the indices without primes refer 
to the two particles of the first Werner pair (the only ones that are 
not tested for spin up) and the primed indices refer to all the other 
particles (that are tested for spin up). The Vuu'u",nn'n" matrix elements 
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of Bob's unitary transformation satisfy a relationship similar to (19). 
The generalization to a larger number of Werner pairs is obvious. 

After the execution of the unitary transformation ([T8|), Alice and 
Bob have to test that all the particles, except those labelled by the 
first (unprimed) indices, have their spin up. They discard any set of n 
Werner pairs where that test fails, even once. The density matrix for 
the remaining "successful" cases is thus obtained by retaining, on the 
right hand side of Eq. (|l5), only the terms whose primed components 
are zeros, and then renormalizing the resulting matrix to unit trace. 
This means that only two of the 2*^ rows of the U matrix, namely those 
with indices 000. . . and 100. . . , are relevant (and likewise for the V 
matrix). The elimination of all the other rows greatly simplifies the 
problem of optimizing these matrices. We shall thus write, for brevity, 

where /i = 0, 1. Then, on the left hand side of Eq. ([T9|) , we effectively 
have two unknown row vectors, Uq and Ui, each one with 2" components 
(labelled by Latin indices mm'm"). These vectors have unit norm and 
are mutually orthogonal. Likewise, Bob has two vectors, Vq and Vi. The 
problem is to optimize these four vectors so as to make the expectation 



value of the Bell operator ||T7 



C ■.= AB + AB' + A'B - A'B', (21) 

as large as possible. 

The optimization proceeds as follows. The new density matrix, for 
the pairs of spin-i particles that were not tested by Alice and Bob for 
spin up (that is, for the first pair in each set of n pairs), is 

N U ^^rnm'm" V^,nn'n" Pmn,st Pm'n',s't' Pm"n",s"t" U^ggigii V^^j/j//, (22) 
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where is a normalization constant, needed to obtain unit trace (A^~^ 
is the probabihty that all the "spin up" tests were successful). We then 
have for fixed pnew and all possible choices of C, 

max [Tr (Cp„ew)] = ^Vm, (23) 

where M is the sum of the two largest eigenvalues of the real symmetric 
matrix T^T, defined by 

Tpg ■= Tr [((Tp (g) O-g) Pnew]- (24) 

(In the last equation, ap and cr^ are the Pauli spin matrices.) Our 
problem is to find the vectors and that maximize M. 

At this point, some simplifying assumptions are helpful. Since all 
matrix elements Pmn,st are real, we can restrict the search to vectors 
f/^ and Vu that have only real components. Furthermore, the situa- 
tions seen by Alice and Bob are completely symmetric, except for the 
opposite signs in the standard expression for the singlet state: 

These signs can be made to become the same by redefining the basis, 
for example by representing the "down" state of Bob's particle by the 
symbol (^i), without changing the basis used for Alice's particle. This 
unilateral change of basis is equivalent a substitution 

yu,nn'n" —>■ { — Vu,nn'n", (26) 

on Bob's side. The minus signs in Eq. (^ also disappear, and there 
is complete symmetry for the two observers. It is then plausible that, 
with the new basis, the optimal Ui, and are the same. Therefore, 
when we return to the original basis and notations, they satisfy 

K,nn'n" = (— l)'^'*'"^" Uu,nn'n"- (27) 
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We shall henceforth restrict our search to pairs of vectors that satisfy 
this relation. 

After all the above simplifications, the problem that has to be 
solved is the following: find two mutually orthogonal unit vectors, Uq 
and Ui, each one with 2" real components, that maximize the value of 
M{U) defined by Eqs. (|23| ) and (|2^). This is a standard optimization 
problem which can be solved numerically. Since the function M{U) is 
bounded, it has at least one maximum. It may, however, have more 
than one: there may be several distinct local maxima with different 
values. A numerical search leads to one of these maxima, but not nec- 
essarily to the largest one. The outcome may depend on the initial 
point of the search. It is therefore imperative to start from numerous 
randomly chosen points in order to ascertain, with reasonable confi- 
dence, that the largest maximum has indeed been found. 



4. NUMERICAL RESULTS 

In all the cases that were examined, M{U) turned out to have a 
local maximum for the following simple choice: 

f/o,oo... = t^i,n... = 1, (28) 

and all the other components of Uq and Ui vanish. Recall that the 
"vectors" Uq and Ui actually are two rows, f/ooo... and f/ioo...; of a uni- 
tary matrix of order 2" (the other rows are irrelevant because of the 
elimination of all the experiments in which a particle failed the spin-up 
test). In the case n = 2, one of the unitary matrices having the prop- 
erty (|28| ) is a simple permutation matrix that can be implemented by a 
"controUed-NOT" quantum gate ]18[ . The corresponding Boolean oper- 
ation is known as XOR (exclusive or). For larger values of n, matrices 
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that satisfy Eq. ( p8D will also be called xoR-transformations. 

It was found, by numerical calculations, that xoR-transformations 
always are the optimal ones for n = 2. They are also optimal for 
n = 3 when the singlet fraction x is less than 0.57, and for n = A when 
X < 0.52. For larger values of x, more complicated forms of Uq and Ui 
give better results. The existence of two different sets of maxima may 
be seen in Fig. 1: there are discontinuities in the slopes of the graphs 
for = 3 and 4, that occur at the values of x where the largest value 
of (C) jumps from one local maximum to another one. 

For n = 5, a complete determination of Uq and Ui requires the opti- 
mization of 64 parameters subject to 3 constraints, more than my work- 
station could handle. I therefore considered only xoR-transformations, 
which are likely to be optimal for x<0.5. In particular, for x = 0.5 
(the value that was used in Werner's original work |^), the result is 
(C) = 2.0087, and the CHSH inequality is violated. This violation 
occurs in spite of the existence of an explicit LHV model that gives 
correct results if the Werner pairs are tested one by one. 

These results prompt a new question: can we get stronger insepa- 
rability criteria by considering p (g) p, or higher tensor products? It is 
easily seen that no further progress can be achieved in this way. If p 
is separable as in Eq. (0), so is p (S> p. Moreover, the partly transposed 
matrix corresponding to p ® p simply is o" Cg> cr, so that if no eigenvalue 
of cr is negative, then a (g) a too has no negative eigenvalue. 
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Caption of figure 

FIG. 1. Maximal expectation value of the Bell operator, versus the 
singlet fraction in the Werner state, for collective tests performed on 
several Werner pairs (from bottom to top of the figure, 1, 2, 3, and 4 
pairs, respectively). The CHSH inequality is violated when (C) > 2. 
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